Pairing imbalance in BCS-BEC crossover of inhomogeneous three-component 

Fermi-gas in two dimensions 
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We in this paper investigate the phase diagram associated with the BCS-BEC crossover of a 
three-component ultracold superfluid-Fermi-gas of different chemical-potentials and equal masses in 
two dimensions. The gap order parameter and number densities are found analytically by using the 
functional path-integral method. The balance of paring will be broken in the free space due to the 
unequal chemical-potentials. We obtain the same particle number-density and condensed fraction 
in the BCS superfluid phase as that in a recent paper (Phys. Rev. A 83, 033630 ), while the Sarma 
phase of coexistence of normal and superfluid Fermi gases is the characteristics of inhomogeneous 
system. The minimum ratio of BCS superfluid phase becomes 1/3 in the BCS limit corresponding 
to the zero-ratio in the two-component system in which the critical point of phase separation is 
esjep ~ 2 but becomes 3 in the three-component case. 

PACS numbers: 03.75.Hh, 03.75.Ss, 05.30.Fk, 74.20.Fg 



I. INTRODUCTION 



Physics related to the crossover from Bardeen-Cooper- 
SchriefFer (BCS) superfluid of Fermi-pairs to the Bose- 
Einstein condensate (BEC) of molecular dimers in ul- 
tracold Fermi gas (FG) has become a central topic in 
both experimental[l|[3|0|[3l and theoretical Q communi- 
ties. The Feshbach resonance plays an important role in 
practical experiments to manipulate the interaction be- 
tween particles in a wide range, which is crucial to form 
the Fermi-pairs and induce the crossover. 

Researches on the population balance and imbal- 
ance in the crossover for nonuniform (the unequal 
chemical potentials or masses) two-component FG in 
2-dimcnsion (2D) [1] d [H and 3-dimension (3D) 
fll][12] 13][14][l^, have attracted considerable atten- 
tions in recent years since unequal particle-number in 
different components may lead to the mismatch of Fermi 
surfaces, where the different components refer usually to 
atoms of different internal states or different isotopes. 
It has been demonstrated that there exist various phe- 
nomena, such as the Sarma superfluidity [l^ with un- 
matched Fermi surfaces and phase separation between 
paired-superfluid and excess-normal fermions depending 
on the population imbalance in the 3D resulted from 
the unequal chemical-potentials and masses. The Sarma 
phase stands for a shell structure of the FG in which 
the core region is superfluid, while the outer region is 
normal [l7|. The phase separation refers generally to the 
spatial separation of unpaired fermions from a BCS su- 
perfluid of equal densities because of the existence of 
pairing gap in the superfluid region(l8|. Thus the Sarma 
phase is also a type of phase separation. 
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2D FG possess a great advantage that the variation of 
bound-state energy in 2D is more suitable to display the 
BCS-BEC crossover than the interaction strength in the 
3D, where a critical coup ling strength is required to form 
a s-wave bound statefT^H^I- Both the energy-gap order- 
parameter, condensate density and population imbalance 
can be expressed as a function of binding energy and the 
phase boundary between normal and Sarma phases is 
identified in the BCS-BEC crossover^ 

The three-component FG have been investigated in 
2D[22] 23][M] and 3D [25] 26] [27] 3] free-spaces respec- 
tively showing the asymmetry of phase diagram by the 
unequal interaction lengths 28]. However, the asymme- 
try is not only induced by the unequal physical quanti- 
ties but also by the pairing itself between different com- 
ponents even in a homogeneous system, which leads to 
population imbalance. It is demonstrated that for a ho- 
mogeneous three-component ultracold FG with a U (3)- 
invariant attractive interaction, the system exhibits pop- 
ulation imbalance in the crossover, which connects the 
BCS regime of pairing state and the BEC limit with 
three species of molecules, since the paring breaks U{3) 
symmetry. The spontaneous population imbalance, first 
noted in Ref. [20] is an important feature of the three- 
component system. Moreover it is shown in Refs. [13] 
and 31] that BCS superfluidity and population imbal- 
ance can coexist in weak-coupling regime of multicompo- 
nent Fermi systems. The BCS-BEC crossover of uniform 
three-component FG with U{3) symmetry was investi- 
gated rece ntlyfaH demonstrating the asymmetry induced 
by paring. The three-component superfluidity has been 
realized already by the lowest three hyperfine spin-states 
in the mixture of ^Li atoms and, moreover, the exper- 
imental analysis of the condensate fraction of ultracold 
two-component FG in the BCS-BEC crossover is shown 
to be in good agreement with the mean-field theoretical 
predictions [3^ [33| and Monte Carlo simulations 0j. 

In a most recent paper (2^ [2^ the analytical formu- 
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las of number densities and condensate fraction for ho- 
mogeneous three-component FG of both 3D and 2D are 
obtained in the BCS-BEC crossover in which, although 
the third component is assumed not to have interaction 
with the other two components the asymmetric interac- 
tions affect the formation of condensation and lead to the 
population imbalance [23| [23| . 

We in this paper investigate the BCS-BEC crossover 
of the three-component FG with asymmetric interactions 



the same as in Ref. [23| and Ref. [2J], but unequal 
chemical potentials in 2D free-space. We are partic- 
ularly interested in the Sarma superfluidity and phase 
separation in the crossover induced by the asymmetry 
of two interacting components. Analytical results of the 
particle-number density and the gap order parameter are 
obtained with the functional path-integral method. Then 
we analyze the condensate fraction based on the phase 
diagram and demonstrate that the Sarma phase is the 
characteristics of the asymmetric system. 



II. PARTICLE-NUMBER AND GAP 
EQUATIONS 

The non-uniform three-component FG, which we are 
going to consider, can be described by the Hamiltonian 

k,i=R,G,B 

+ (5fl,G^I+q/2,ijV^!.k+q/2,G^-k/-Hq/2,GV'k/+q/2,fl 
k,k/,q 

+ 5S,flV'I+q/2,fl^Ik+q/2.S^-k/+q/2,BV'k/+q/2.fl, 
+ 5S,GV4+q/2,BV'-k+q/2.GV^-k/ +q/2,GV'k/+q/2,B) 

(1) 

where '0^ ^ ('0k, i) is the Fermi-field creation (annihilation) 
operator creating ( annihilating) a particle of color-index 
i — R,G,B (denoting red, green and blue respectively) 
and momentum k. ^k,i = — with Ek — k"^ I (2m) 
being the kinetic energy. The three components possess 
the same mass m but unequal chemical potentials fii. 
g with a negative value is used to describe the attrac- 
tive interatomic interaction. The rest three sums in the 
Hamiltonian denote the pairing among particles of dif- 
ferent colors while there is no pairing in the same com- 
ponent of FG. We furthermore assume that the pairing 
takes place only between two components, say, the red 
and green particles. Thus the Hamiltonian is reduced to 
a simple form 



H = Y Ck,i'0k,»V'k,i 
k,i=R,G,B 

+ J2 5i?.GV'l+q/2,iiV'Ik+q/2,G'/'-k/+q/2,GV'k/+q/2,fl 



k,k/. 



By using the functional path- integral treatment, we 
obtain the thermodynamical potential of saddle-point, 
which can be written as[35l| 



A2 

+ £,KG-E^,.- + fF{^iB)^KB ] (3) 

9 

where ii^k,± = (Ck + + A^)^/^ ± ^k,- are the quasiparticle 
and the quasihole energies with ^k,± = ('Ck,fl±Ck,G)/2 = 
k'^ /{2m) — ^± and energy-gap A is the order-parameter 
of FG. The reexpressed chemical potentials ^± = {^ib. ± 
Plg)/2 indicate the chemical potential imbalance between 
red and green particles. The particle number of color- 
index i is obtained from the saddle point approximation 
as A^o.i = —^^o{^,^J^R,^J'G,^J^B) /d^i, and the explicit 
formulas of the three-color particle-number densities rii = 
^o,i/^ (* = G, B) are given respectively by 



riR 



riG 



where 



/ -^[i^-fiEK-))vl + fiEK+>l] 
J (ztt) 



(2^) 
(2^) 



(2n 



2 [(l-/(i?k.+))«: 

2 [fP ifiB)] 



HEk- 



(4) 

(5) 
(6) 



;i-f a,+/£^k,+)/2, 



(1 — ^k,+/£'k,+) /2 and the Fermi thermal distribution- 
function has the usual form f{x) = (exp(/3a;) -1-1) ^ . 
The number equations ^ and ([5]) of paired compo- 
nents must be coupled with the gap equation given by 
SSo/SA = 0, ( S'o = pno), which is 



f{Eu 



f{EK-) 



(7) 



k,+ 



The binding energy e^, which is more suitable than in- 
teraction strength g to describe BCS-BEC crossover as 
demonstrated in the introduction, can be determined 
from the bound-state energy equation 



- = -y 



1 



(8) 



(2) 



The particle-number (|4]), ([5]), dUand gap equation ([7]) 
are going to be solved self-consistently under the zero- 
temperature limit that 6'(— i?k,±) = lim^_^oo /(£'k,±), 
where 9{x) is the Heaviside function and the analytic 
results are shown in the following section. 



III. CRITERION OF THE PHASE TRANSITION 

At zero-temperature the saddle-point self-consistent 
equations dH), ©, (H]) and (O are suitable to describe 
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the BCS-BEC crossover, from which we can calculate the 
particle number and the binding energy e^. When the 
quasiparticle energy vanishes -Ek;,+ = 0, the zero-point 
energies e± are found as 



1 

2 

^1+ ± J fi 



4A2 



(9) 



where we have used the property of Fermi thermal 
distribution- function that f{x) — at x < and f{x) = 
1 at X > at zero-temperature with the integral vari- 
able, i.e. the kinetic energy — fc^/ (2m). In the chem- 
ical potential range — A < /i_ < A, the zero-point ener- 
gies e± are complex values, so the quasiparticle energies 
i?k,+ and £^k,- are positive indicating the BCS super- 
fluid state. When > A, e± are real and moreover 
in the branch of e± > 0, the quasienergy i?k.+ is nega- 
tive, while the kinetic energy in the region e_ < et < £+ 
the quasienergy i?k,- is always positive. To have the 
kinetic energy Sk being positive, its value-range should 
be restricted to [0,e+] when e_ < giving rise to the 
Sarma phase, which is the coexistence of the normal and 
superfluid FG. 

The particle-number of EEC phase is strictly re- 
lated to the off-diagonal long-range order (ODLRO) 
in interacting many-body systems, which was pro- 
posed in terms of the asymptotic behavior of two- 
body density matrixfsll defined as p(r^, Tjj, r^j, rg) = 

(^fi('"fl)^G(rG)V'fl(ri?JV'G(rG)) , where ip,{r) and V'l(r) 
are the fermionic creation and annihilation operators of 
i-th species in the position r respectively and (...) stands 
for quantum-mechanical expectation value. For a homo- 
geneous FG, the ODLRO of Fermi superfluid is char- 
acterized by the asymptotic behavior of density matrix 



Pi^R^^G^^R^^c) for 



re 



oo. While 



the largest eigenvalue Nc of the two-body density matrix 
gives rise to the number of Fermi pairs in EEC phase, 
which is denoted by 

N, = j d^VRd^rG\{i^R{vR)^G[-rG))\^ 
and it is straightforward to showfs?'] that 



N, 



(2^) 



2 2 



(10) 



In this paper, we moreover define a physical quantity 
called ratio of unpaired particle-number imbalance as: 



UR - ng + nB 
nR + UG + riB 



(11) 



which with value in the range [0, 1] can identify the phase 
transition. 



A. Sarma superfluid phase 

At first, we consider the zero-point energy range 
[£_,£+], which can be classified into three cases that 
< e_ < e+, £_ < < s+, and £_<£+< 0. The 
third case, which is the same as in the chemical potential 
region —A < /i_ < A, is the ECS phase. The two other 
cases are both the Sarma phase. When -Bk,+ < with 
k ^ and < £_ < £+, the Sarma phase is in the ECS 
regime while in the EEC regime for fc = 0, £_ < < £+. 
The particle-number density equation is 



n = riR + ug + nB — mep/Tr 



(12) 



where ep is the Fermi energy. The number-density dif- 
ference between red and green particles of the Eqs. ([4]) 
and ([5]) is seen to be 



ur-ug = m/ (27r) (£+ - £_) 



(13) 



Thus particle- number densities in the range [£_,£+] are 
found from Eqs. (|4|) and (O as 



riG = 

nR = m/ (27r) (£+ - £_) 
ub = m/is/ (27r) 6 {(Ib) 

along with vanishing gap A = 0. For the sake of simplic- 
ity, we let the chemical potential of blue component be 
fiB = IJ'+ throughout the paper. 



1. Sarma phase P2 in the BCS regime 

For the case one (0 < £_ < £+) with i?k,+ < , i?k,- > 
0, we have fp {Ek.,+) — ^ , fr (Ew,-) = the particle 
number density of green component evaluated from the 
Eq. (m is given by : 



HG 



^(1/2- . Jr':l.., ) (14) 



2V(£-M+)'+A2 

With the result of Eq. ^ together we have 

nR-nG = (15) 

The energy-gap and bound-state energy are determined 
by the equations 




mds 



1 



An 



^(£-^+)VA2 



1 
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ide 



4tt e + eb 



(16) 



(17) 
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which can be combined together as 

\ mde 1 







1 mde 1 



2 7o An e + EB 



(18) 



where cj± = ^,2^0 = ^ ; 2;^ = ^ . Working out the 



integration in Eq. (|14p . the particle number density has 
the form: 



mA . f. ,2 -, 

nG=— [w- - (w+-a;o) +1 



71_B = (27r) (^+) 



(19) 
(20) 



and the gap parameter is obtained from Eq. (|18l) 



A 



it can be easily rewritten as 



A_ 







(uj+-xo)+\J {uj+-xof+l 









(w+-.To)Vl- \J {uj--xof+l+XQ 9{xa )+xo+^/ x^+l 

(21) 

For n = nil + + "-s = mep/T^ with kp — \J~^^, the 
binding energy eg is seen to be 




and the ratio of unpaired particle-number imbalance P 
is easily obtained from Eq. ([TT|) as: 

a;+^_+a;o6'(a;o) 



(a;+-a;o)Vl- (w.-xq) Vi+xq 6^ (xq )+xo+ VaJg+T 

(23) 

The condensate fraction nc= Nc/V is found from the Eq. 
(fTOl) as 



Am 
'"8^ 



^+arctan (^)-2 arctan { -I 



(24) 



^. Sarma phase PI m BEC regime 

In case two (e_ < < £+) the integral range is from 
to e+ where the minimum value of £'k.+ is located at 

= 0, so the system is in the BEC regime. The number 
density and energy-gap can be shown from the Eqs. (|13p, 
(P, and ®as: 



Ztt 



no 



mA 

An 



\J [uj+ - xof + 1 - a;+ + xq 



(25) 



(26) 



A 



\J {uj+ - xq)^ + 1 + a;o6'(a;o) + a;o 
The binding energy is 

(w+ - xq) + (w+ - xaf + 1 



\J (iu+ - xo)^ + 1 + .xo6'(.xo) + xo 



and 



+ xo6'(xo) 



{io+ - xo)^ -I- 1 -I- Xo6'(xo) -f- xo 
We can easily get the condensate fraction 



Am 
"8^ 



arctan 

2 



^1 



(27) 



(28) 



(29) 



(30) 



B. BCS superfluid phase 

For positive -Ek.it and £_<£+< 0, particle number 
densities of the red and green components are equal seen 
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FIG. 1: (Colour online) The number-density ratios nn/n, and 
nB/n, as a function of the bound-state energy es/^F- 



from the Eqs. ([4]), ([5]) and the chemical potential is in 
the region —A < < A. Number densities and gap 
parameter are given by: 



ur = no 



A 



mA 

47r 



x2 + l 



(31) 



(32) 



a;o6'(a;o) + a^o + ^/xl + l 
In this case, the binding energy can be shown as Eq. (|22 



and 



P 



1 - 



a;o6'(a;o) 



a;o6'(xo) + + V^^ + l 



(33) 



(34) 



The ratios of particle nmxiber densities obtained from 
Eqs. (|3T|) and ([20l) are given by 



n 



n 



1 Xq + yjx\ -h 1 



2 a:o6'(xo) + 2:0 + V^zif + T 



and 



n 



xqQ{xq) + + yirfTT 



(35) 



(36) 



which in agreement with the Ref. (23| are displayed 
in Fig.l, where we see the balance of number densities 
nR/n = nc/n = ns/n = 1/3 .in the BCS regime. With 
increase of bound-state energy cb/^f, the number den- 
sity of red particles UR/n increases while nsln decreases. 
At the energy value cb/^f =3, the red density-ratio be- 
comes 1/2 and the blue one is zero. The red and green 
density-ratios are able to sustain 1/2 with higher value of 



eb/cf in the BEC regime. The condensate fraction Eq. 
([TU]) gives rise to the expression 



n 



1 



-I- arctan (xq)] 



: a;o6'(a;o) + a^o + ^/xl + 1 



(37) 



Phase transition and boundary 



The frequency uj± — xq = ± i/j/g — 1 with yo — fi^/A, 
is a characteristic quantity of the system at hand with 
which the ratio P and bound state energy cb/^f can be 
expressed as the function of two parameters Xq and j/o- 

In the Sarma phase ratio P and bound state energy 
^bI^f become functions of parameters xq and yo 



P 



xoe{xo)+2^yl - 1 



xo9{xo) +Xo + yjxl + l 



EF 



3 


^Jxl + I-X0 




\Jvl - 1 + yo 




xoO{xo) + xo + ^xl + l 




yo - V^o - 1 



(38) 



(39) 



in the region < e„ < £+, and 



P 



xoOjxo) +XQ + y/y^-1 
xoO{xo) + a;o + yo 



CB _ 3 [Vyo - 1 + yo 

e_F XQd{xo) +xo + yo 



(40) 



(41) 



in the case of e„ < < e+ respectively. The phase 
transition and boundary are determined by the number 
of zeros of quasienergies i?k,± in the momentum space. 
Sarma phase in two ranges with one and two zeros respec- 
tively, which identify the normal gas, does not exist in the 
three-component FG of equal chemical potentials. The 
BCS superfluid phase in two ranges has no zero, while 
two negative zero-point energies correspond to negative 
quasienergy also different from the standard BCS. 

In the phase diagram of Fig. 2 obtained from Equations 
of P and sb/^f, -PI stands for the Sarma phase with 
effective Fermi-surfaces formed by the two zero-points 
of excitation-spectrum £'k.+ and one negative zero-point 
energy in the condition — ^?/q — 1 < xo < \/y'o ~ 1- 

The Sarma phase P2 has effective Fermi-surfaces re- 
sulted from two zero-points of £'k.+ -in the condition 
2^0 > \/y§ — 1- The BCS phase indicated in Fig 2 has 
no effective Fermi-surface of zero quasienergies. More 
specifically in the condition yo < 1- -E'k.+ has no zero- 
point in one case, which is consistent with the classical 
BCS superconductivity, and two negative zero-point en- 
ergies in the oth er case with the condition xo < and 
Xo < -\/vl " 1- 
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suppress the formation of long-range phase coherence at 
non-zero temperature |38i] and remain non-neghgible even 
at zero temperature. It has been demonstrated the quan- 
tum fluctuations kee p on playing a crucial role also when 
imbalance is present [39|. As a consequence, in order to 
describe Fermi gases in the whole range of the BCS-BEC 
crossover at zero temperature the mean-field approxima- 
tion might be not sufficient, and a more accurate account 
of quantum fluctuations is necessary. The research be- 
yond the mean-field approximation is under development 
referring to the BCS-BEC crossover. 

IV. CONCLUSIONS 



FIG. 2: (Colour online) Phase diagram in the plane of en/cF 
and P. The Sarma phase is divided to two parts PI and P2 
by the critical line P = 3 — es/eF- 

From Fig. 2, it is found that at the BCS limit (es/eF = 
0) P = 1/3 in both the BCS and the Sarma phases, which 
is consistent with the Fig.l at es/e^ = and is different 
from the two-component FG (in which P = at the 
BCS limit) due to the existence of the third component. 
The critical point a;o = 0, where es/^F = 3 and P — 
0, is a common point of the three phases. The Sarma 
phase is divided to two parts PI and P2 by the critical 
line P = 3 — es/ei?, while the two parts are divided by 
the line ^bI^f = 2 in two-component FG. In the BCS 
superfluid phase P = for es/^F > 3 and it decreases 
from P = 1/3 to P = with the increase of es/ei? in 
the region tsl ^f < 3 different from the two-component 
FG in which the BCS phase goes along P = in the 
P — eb/^f plane. Nl and N2 in Fig. 2 respectively denote 
the partially and fully polarized normal phases, which 
are also abundant phases and do not exist in the three- 
component FG of equal chemical-potentials. With the 
increase of P and es/eF, the phase transitions of BCS, 
normal, Sarma phases P2, PI take place in order. When 
es/cF > 3, there are only red and green components in 
the crossover and the BCS superfluid phase emerges at 
P = i.e. the balance pairing nji = no, in consistence 
with the Fig.l, which shows nc/n — riR/n = 1/2. 

For a Fermi gas in 2D, the quantum fluctuations play 
a more crucial role with respect to the 3D case, which 



The pairing physics in BCS-BEC crossover is stud- 
ied with particular interest in the chemical-potential 
asymmetry induced imbalance. The analytical energy 
gap, particle-number density and the condensate fraction 
of ultracold three-component FG of unequal chemical- 
potentials in 2D are shown in exact agreement with the 
results in Ref. [1^ of homogeneous system, while Sarma 
phase, which is divided into two parts P2, and PI, par- 
tially and fully polarized normal phases are new and are 
resulted from the asymmetric chemical potentials. Phase 
diagram and transitions in the P — cb/^f plane are also 
compared with the inhomogeneous two-component FG. 
Although the third component is assumed not to pairing 
with the other two components the ratio P, which van- 
ishes in the BCS limit for the two-component FG[9], be- 
comes 1/3 in the three-component system with the crit- 
ical point of phase separation eb/^f = 3. 
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